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ABSTRACT 



Context. Previous studies have shown that many post-AGB stars with dusty disks are associated with single-lined binary stars. The 
inferred orbital separations are too small to accommodate a fully-grown AGB star, therefore, these systems represent a new evolu- 
tionary channel that bypasses a full AGB evolution. 

Aims. To verify the binarity hypothesis on a larger sample, to reveal the nature of the companions, and to probe the disk structure and 
eventually the disk formation mechanisms, we started a high-resolution spectral monitoring of ~40 field giants whose binarity was 
suspected based on either a light curve, an infrared excess, or a peculiar chemical composition. 

Methods. Starting from the spring of 2009, we monitor the program stars with a new fiber echelle spectrometer HERMES. We mea- 
sure their radial velocities (RVs) with a precision of ~0.2 km s"', perform detailed photospheric abundance analyses, and analyze the 
time-resolved high-resolution spectra to search for line-profile variability. 

Results. Here we report on the discovery of the periodic RV variations in BD+46°442, a high-latitude F giant with a disk. We interpret 
the variations due to the motion around a faint companion, and deduce the following orbital parameters: P^^ = 140.77 ± 0.02'', e = 
0.083 ± 0.002, a sin ; = 0.31 AU. We find it to be a moderately metal-poor star ([M/H] = -0.7) without a strong depletion pattern in 
the photospheric abundances. Interestingly, many lines indeed show periodic changes with the orbital phase: Ha switches between a 
double-peak emission and a P Cyg-like profiles, while strong metal lines appear split during the maximum redshift. Similar effects 
are likely visible in the spectra of other post-AGB binaries, but their regularity is not always realized due to sporadic observations. 
We propose that these features result from an ongoing mass transfer from the evolved giant to the companion. In particular, the blue- 
shifted absorption in Ha, which occurs only at superior conjunction, may result from a jet originating in the accretion disk around the 
companion and seen in absorption towards the luminous primary. 



Key words. Stars: abundances - binaries: spectroscopic - Stars: individual: BD +46° 442 - 
and post-AGB - ISM: jets and outflows 
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1. Introduction 



Recentlv lvan Aarle et al.l (1201 ih investigated the spectral energy 
distribution (SED) of optically bright candidate post-asymptotic 
giant branch (post-AGB) stars in the Large Magellanic Cloud. 
Surprisingly, they found that in about half of candidates, the SED 
indicates that the circumstellar material, which is the relic of the 
AGB mass loss, is not in the form an exanding envelope, but 
rather in the form of a stable, likely Keplerian, disk. 

In our ow n Galaxy, post-A GB stars with disks are common 
as well (e.g. Ide Ruvter et al.l 12006). but statistics of Galactic 
post-AGB stars are less complete, and the lack of well con- 
strained distances (hence, luminosity) makes the interpretation 
of these numbers more difficult. Galactic sources, however, have 
an obvious advantage of being closer and easier to study. 

The leading hypothesis explaining the envelope bifurcation 
into the outflows and the disk types, is that disks form only 
around binaries. The observational foundation comes from the 



* Based on observations made with the Mercator Telescope, operated 
on the island of La Palma by the Flemish Community, at the Spanish 
Observatorio del Roque de los Muchachos of the Institute de Astrofisica 
de Canarias. 



fact that indeed, many disc s ources in our Galaxy turn ed out to 
be single-lined binaries (e.g. Ivan Winckel etalll2009l and ref- 
erences therein). Presumably disk formation happens through 
an interaction with the companion when the more massive 
component becomes a red giant (RG) or an asymptotic giant 
branch (AGB) star (Ivan Winckell 20031). The process itself is 
poorly understood. iLivio & Sokei ( 19881 ) and ISandquist et al.l 
(1998), among others, showed that a circumbinary disk could 
be a remnant of a common envelope created when the Roche 
lobe of an AGB star engulfs a companion. Other s cenarios 
were proposed by e.g. iMastrodemos & MorrisI (Il999l) . where 
a disk-like structure is created in a wind-accretion scenario or 
perhaps as a result of the AGB wind and the accretion jets 
interaction(Ak ashi & So ker 2008). 

The impact of the disk onto the system cannot be underes- 
timated. Many of the Galactic and even LMC post-AGB disk 
objects show a peculiar chemical photospheric composition that 
resembles the gas phase of the interstellar medium (ISM): el- 
ements of higher condensation temperature, such as ir on, are 
less abund ant than volatiles (e.g. Giridh ar et al . 2005; Ma as et al.l 
2005; Gi elen et aH2 009; Sumangala Rao et al.l2012l) . This is be- 
lieved to be a result of pollution from a disk, which may have 
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already dissipated in some systems and therefore has no foot- 
print in the SED. In the disk, refractory elements nucleate into 
dust grains, which are segreg ated from the ^as, and the depleted 
gas falls back onto the star (IWaters et al.lll992 ). On one hand, 
this effect helps to identify current or former disk systems, but 
on the other hand masks the original metallicity of the stars and 
the eventual products of the AGB dredge-ups, complicating the 
identification of the evolutionary stage. 

At a later evol ution ary stage, as was proposed by 
IWaelkens et alj d 19961) and iJorissenI (Il999l) . a disk system may 
turn into a Barium star, when the current post-AGB primary 
evolves into a white dwarf (WD) and supposedly a main se- 
quence (MS) companion into a red giant. Barium stars are bi- 
naries where a red giant bears signs of contamination by the 
nuclear-synthetic products (such as Barium) by the former AGB 
comp anion (now a WD). By adoptin g the form alism of Lubow 
(1201 0. and references therein), Der mine et al.l (tZOlO ) explored 
the binary-disk interaction regime relevant in Barium stars and 
showed that circumbinary disks may be important at pumping 
the eccentricities of Barium stars to the observed values. 

The Galactic objects with known orbits cannot have evolved 
on single-star evolutionary tracks: typical AGB star has a radius 
~ 1 AU, while many post-AGB binaries with a disk have smaller 
separations, which brings a possibility that these systems started 
interacting already at the red giant stage. Some of the objects 
may have never gone through the AGB stage, unless the inter- 
action resulted in the orbit shrinkage. With the observed peri- 
ods, the full AGB evolution must be shortcut and this is cor- 
roborated by the finding that the majority of disk sources have 
oxygen- rich photo spheres a nd dis ks dominated by crystalline 
silicates dGielen et al. 20081 l201lh . Another ambiguity pertain- 
ing the evolution of these systems is whether the stars will still 
show up as planetary nebulae (PNe), with the circumstellar ma- 
terial locked in a disk. While binarity helps to explain bipolar 
PNe, the periods of binary stars in PNe are at least tw o orders of 
magn itude shorter than in the disk post-AGB systems (iDe Marcol 
I2009h . 

The theories of disk formation cannot predict timescales 
for disk re-accretion and dissipation either The insight can be 
gained by gathering information on the inner-most regions of 
the disk from optical and near-IR spectroscopy. Previous spec- 
troscopic studies of post-AGB stars explored mainly stellar ra- 
dial velocities (RVs) and the basic properties of both the pho- 
tosphere and the disk. Only briefly did they touch on the signs 
of activity at the star-to-disk interface, such as emission in Ho- 
or circumstellar absorptions in the Ca II H&K and Na I D lines, 
and the interpretation was usually that these features trace a 
subsiding AGB wind. The notion was that after a circumbinary 
disk is formed, the system enters a passive evolution as a de- 
tached binary. No success has been achieved at uncovering gas 
re-acc retion (with only on e tentative detection of the winnowing 
gas bv lHinkleet al]|200l . In none of the systems, the compan- 
ion is detected so far. To explore the details of the disk-star in- 
teraction, systematic time-resolved investigation of the disk and 
chemically-peculiar post-AGB systems is therefore highly desir- 
able. 

A literature compilati on of Galactic p ost-AGB s was assem- 
bled in the Toruh catalog (ISzczerba et al.l l2007fl and consists of 
~400 objects. I n this catalog, 5 1 sour ce have been identified as 
disk sources in Ide Ruvter et al.l (l2006h . and this number has in- 
creased ever since. About two dozen post-AGBs of De Ruyter 
et al's list have been identified to be single-lined spectroscopic 
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binaries with periods from 115 to 2600 days and orbital sepa- 
rations from 0.1 to 3.6 AU. The latest addition can be found in 
van Winckel et al. (2009) who photometrically and spectroscop- 
ically monitored a selection of low-amplitude pulsators from the 
southern hemisphere. All low-amplitude pulsators are binaries, 
but the suspected binary nature of the remaining 60% of De 
Ruyter et al's list remains to be established. This is an obser- 
vational challenge, because many objects show large amplitude 
pulsations making the in terpretation of rad ial velocity variations 
not straightforward (e.g. lMaas et al 1 l2002h . 

In the spring of 2009, we started a high-resolution spectro- 
scopic monitoring programme on 34 Galactic post-AGB can- 
didates with disks and a few objects without the IR excess, 
that are observable from the northern hemispher e, using a 
newly built optical echelle spectrograph HERMES ( Raskin et al.l 
201 1). These objects were required to satisfy at least two of the 
following criteria: an IR excess characteristic of a disk, an RV 
Tau-like variability, variability suggesting a possible obscuration 
by a disk, low surface gravity combined with the location away 
from the Galactic plane, deficiency in the refractory elements. In 
this paper we report on one of the first new short-period binaries 
that we uncovered after three years of observations. BD-i-46°442 
(SAO 37487, IRAS 01427-H4633) is a poorly-studied V = 9.5 
high-latitude Galactic F gianj^ {b - 15°), not known to be vari- 
able in RV or photometry. It was selected in our survey due to 
the strong infrared excess characteristic for disk objects, as orig- 
inally measured by the IRAS satellite. The paper is organized 
as follows: after description of the data (Sect.|2]i, we present re- 
sults of the abundance analyses of BD-i-46°442 (Sect. O, then 
model its SED (Sect.©, the radial velocity (RV) curve (Sect.©, 
and examine the behaviour of the hydrogen and some metal lines 
(Sect. |6]l. In the Discussion section, we compare the discovered 
periodic spectral variations with those in other post-AGB sys- 
tems and in interacting binaries, and discuss theories capable to 
explain this behaviour. 

2. Observations and data reduction 

Each target in our on-going post-AGB survey is observed about 
twice per month during the visibility season. This program is 
part of a much broader effort to study binary interaction pro- 
cesses in evolved stars (Van Winckel et al. 2010). The main goal 
is to measure RVs to detect and characterize the orbital mo- 
tion. For brighter targets, additional studies are performed, such 
as determination of abundances, rotation velocities, and the cir- 
cumstellar features. The spectra are obtained with the fiber-fed 
echelle spectrograph HERMES, attached to the Flemish 1.2 m 
telescope Mercator on La Palma, Canary Islands. The spectro- 
graph is optimized for high resolution, stability, and broad wave- 
length coverage. This is achieved primarily by implementing an 
imag e slicer, an anti-frin ge CCD coating, and a thermal enclo- 
sure (lRaskinetal.ll20Ilh . 

Sixty spectra of BD-i-46°442 were obtained during three sea- 
sons of observations in 2009 - 2012. We used a fiber configu- 
ration that provides a resolution R - 80, 000 - 90, 000 over the 
wavelength range AA = 3800 - 9000A. The exposure times var- 
ied between 500 - 1800 s while the signal-to-noise ratio (S/N) 
between 20 - 90 (as measured near 6500 A). Nightly calibra- 
tions consisted of two sets of biases, flats, and the Th-Ne-Ar 

- When referring to BD+46°442-like systems, the term "giant" is 
used to indicate a low surface-gravity, extended, post-MS object. The 
exact luminosity class and the evolutionary stage can only be estab- 
lished when the parallax information becomes available. 



2 



Gorlova et al: A new evolved binary with a disk 



arc spectra obtained in the evening and in the morning. To re- 
duce a particular object frame, the nearest set was used. At least 
one lAU radial velocity standard was observed nightly and oc- 
casionally a series of arcs, to monitor the fluctuations in the 
zero-point of the wavelength calibration. These tests showed that 
the typical uncertainty on the absolute value of the RV is ~Q.2 
km s"' and is mainly caused by the pressure fluctuations in the 
instrument room (Raskin et al. 2011). The data reduction was 
performed with a dedicated Python-based pipeline, that outputs 
extracted, cosmic -ray cleaned, wavelength calibrated, and order- 
merged spectra. 



3. Atmospheric analysis 

3.1. Atomic data and model atmospheres 

Our atmospheric analysis of BD-i-46°442 is based on the com- 
parison with synthetic profiles of hydrogen lines and on the 
equivalent-width (EW) measurement of metal fines. The lat- 
ter we re obtained using tfie DECH20 code by Galazutdinov 
(ll992lPL tfiat provides a convenient interface for the continuum 
normafization, and the EW measurement by means of a Gaussian 
fit or a direct integration. Line identification was carried out us- 
ing the following on-line resources: Spectroweb stellar atlas by 
^^el (2008, 201 ifl for A < 6800 A, and the Solar atlas by 
iDel bouille et aL, (,1988 r| in combination with the Atomic line 
database by R Van HooQ for the longer wavelengths. 

The line list used for the abundance determination 
was provided by V. Kovtyukh (2010, private communica- 
tion) and is mainly based o n the version published in 
iKovtvukh & Andrievskvl (Il999h . We also follow the proce- 
dure for the determination of the surface gravity (log^) and 
the micro-turbulent velocity (Vtur) from the latter work. Our 
choice is dictated by the fact th at yellow supergiants studied by 
iKovtvukh & Andrievskvl d 19991) spectroscopically are very sim- 
ilar to post-AGB stars. To convert EWs to abundances, we use 
the local thermodynamical equilibrium (LTE) radiative transfer 
code MOOG by C. SnedeiQ with the ATLAS9 model atmo- 
sp heres that include t he up dated opacity distribution functions 
of lCastelfi & Kuruczl(l2003lB 

IKovtvuk h & Andrievskvl ([1999') obtained the bulk of their 
oscillator strengths (loggf-s) by adjusting the laboratory values 
in such a way, that they could reproduce the solar abundances 
of iGrevesse et al.* ('1996) from the EWs measured in the solar 
spectrum of [Kurucz et al. (1984). They used a canonical solar 
model of 5777/4.438/1.0 (Feff/ log;?/ V^J from the set of the 
original ATLAS9 models of Kurucz ( fl992 |) and his WIDTH9 
LTE code 0. We found that this combination of the older at- 
mospheric models and the code tends for the lines of a given 
element to produce systematically slightly different abundances 
than our combination of the atmospheric models and code, both 
for the Sun and BD-i-46°442. Therefore, when giving abun- 
dances of BDh-46°442 relative to the Sun ([X/H]), we introduced 
small corrections (up to +0.08 dex) in the solar abundances of 
[Kovtvukh & Andrievskv ( 1 999) to compensate for this differ- 



http://www.gazinur.com/DECH-software.html 

http://spectra.freeshell.org/spectroweb.html 

http://ljr.bagn.obs-mip.fr/observing/spectrum.html 

http://www.pa.uk;y.edu/~peter/atomic/ 

http://www.as.utexas.edu/~chris/moog.html 

http://wwwuser.oat.ts.astro.it/castelli/grids.html 

http://wwwuser.oat.ts.astro.it/castelli/sources/width9.html 



3.2. Photospheric parameters 

Since there is no information on the reddening or the photo- 
metric variability in BD+46°442, we decided to derive all pho- 
tospheric parameters spectroscopically. For this we analyzed a 
S/N~130 spectrum obtained by averaging the best quality spec- 
tra numbered 3 and 4 from the same night (Fig. [Q. The visual 
examination of all the spectra did not reveal any obvious varia- 
tions in the eff'ective temperature (Tgff). Variations in the cores 
of strong lines have been detected, which we attribute to the cir- 
cumstellar matter (Sect. |6]l, but those lines were not used in the 
EW analysis (we only used lines with EWs < 170 A ). This sin- 
gle spectrum, therefore, should give a good representation of the 
basic photospheric parameters of our star. 

First, we confirmed a spectral type F from Simbad by com- 
paring our spectru m to the spectral libr aries of R. O. GrayH and 
of the VLT-UVES (Bagnul o et ani2003lM Then we computed 
from our EWs the iron abundance, by adopting a range of model 
parameters expected for an F I-III star: T^ff = 6200 - 8000 K, 
log^ = 0.5 - 4.0, and Vu,r - - 12.0 km s"', first using the 
ATLAS9 models with the solar metallicity ([M/H] = 0) and later 
with [M/H] - -0.5 to match the subsolar metallicity obtained in 
the final iteration. To perform an excitation analysis, we plotted 
the iron abundances against the EW and the excitation poten- 
tial of the lower level (x), separately for the Fe I and Fe II lines, 
to examine possible trends. Correct (LTE) atmospheric param- 
eters are obtained when all the lines give the same value of the 
abundance. Trends with the EW can be removed by tuning Vtur, 
while agreement between the average abundances derived from 
the Fe I and Fe II lines (the ionization balance) indicates a correct 
log g. Iron lines are usually used in this method because they are 
the most numerous, but the same values of Vtm and log g should 
in principle hold for all the elements. 

At sub-solar metallicities and gravities, however, neutral 
species become underpopulated compared to the LTE case, due 
to the UV over-ionization and the decreased rate of collisions, 
resulting in the positiv e corrections that need to be added to 
the LTE abundances ( | Rentzsch-Holmlll996t lBergemannll201 lb 
Mashonk ina et alj|201 lb . Observationallv this is manifested in 
the fact that in supergiants, Fe I lines require smaller values of 
Vtur and \ogg to bring their abundances in agreement with the 
better behaved Fe II lines (which is the dominant ionization state 
in F-G stars), as shown for example in IKovtvukh & A ndrievskY 
(1999) for the classical Cepheid 6 Cep and in Ijakeda et aL 
(2007) for some post-AGB stars. We find the same eff'ect for 
BD-i-46°442. Following recommendations in the literature, we 
use the Fell lines for the determination of Vtm, while for the 
determination of log g we extrapolate abundances from the Fe I 
lines to EW = 0, as weak lines are less susceptible to non-LTE ef- 
fects. Following this procedure, we found a correlation between 
the adopted Teff and the derived logg/[Fe/H]. To break this de- 
generacy, we turn to the Balmer and Paschen hydrogen lines. 

Hydrogen lines have the advantage over the metal lines that 
to first order they do not depend on the metallicity. Furthermore, 
we can consider the shape of the extended wings in addition 
to the line depth when comparing to the synthetic profiles. We 
compar ed Balmer lin es in BD+46°442 to the synthetic pro- 
files of IC oelho et al. (2005|3 and Paschen lines to those of 
iMunari & CastellL (.2000.Fl Ha was not considered, as it shows 



http://ned.ipac.caltech.edu/level5/Gray/Gray_contents.html 
' ' http://www.sc.eso.org/santiago/uvespop/ 

'2 http://www.mpa-garching.mpg.de/PUBLICATIONS/DATA/SYNTHSTELL 
http://vizier.u-strasbg.fr/viz-bin/VizieR-4?-source=III/238 
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Fig. 1. Sample of the spectrum of BD+46°442 obtained on 
2009/08/17 compared to the HERMES spectrum of Polaris. 
While Hy6 is similar in both stars, indicating similar tempera- 
ture and gravity, metal lines are noticeably weaker (and broader) 
in BD-l-46°442, indicating a sub-solar metallicity. 



Fig. 2. Hydrogen lines in BD-h46°442 (black dots) compared to 
the best fitted of the available synthetic profiles (grey solid lines), 
that were computed with Ten - 6250 K, log^ = 2.0, [M/H] = 
-0.5. To match the resolution of the model Paschen lines, the 
observed profiles were convolved to R = 20, 000, while B aimer 
lines are shown at the native HERMES resolution of 85,000. 



a strong non-photospheric contribution at all phases. From the 
Paschen lines we chose to consider only Pal4 and Pal7 as the 
least blended ones. We noticed that the Paschen break in the 
models is too deep compared to BD+46°442 and other super- 
giants we have studied; therefore, before making a comparison, 
we had to normalize the model and the observed profiles at the 
pseudo-continuum regions between the neighboring lines. 

In the parameter space considered here, model hydro- 

[ en lines strengthe n both with higher T^s and lower \ogg 
Fremat et al.l Il996t iMunari & Castellil l2000l) . Fortunately, a 
stronger sensitivity of Paschen Unes to gravity allowed to ex- 
clude log^ > 3.0, and > 6500 K along with it. No Paschen 
models are available for log g < 2 for the range of temperatures 
considered here, but they would require Teff < 6000 K, which is 
too low to fit Balmer lines. We therefore adopt T^ff = 6250 + 250 
K, logg < 3.0 (Fig.|2]), and return to the iron lines to put a bet- 
ter constraint on log g. Our estimate is further confirmed by a 
good match of Paschen lines and the shape of H/3 wings between 
BD-i-46°442 and Polaris (Fig. [TJ- Polaris Aa is a well-studied 
short-p eriod Classical Cepheid, for which e.g., lUsenko et alj 
(I2OOI infen-ed Teff = 6000 ± 170 K, logg = 2.2 + 0.3, and 
[Fe/H] = -i-0.07. 

Finally, Tgg is also frequently derived from the requirement 
that the iron abundance derived from the Fel lines does not 
depend on x- Applying this criterion to BD-i-46°442, we ob- 
tain Tgff «i7500 K, which is inconsistent with the hydrogen fits. 
However, as we discussed above, Fe I lines may not be reliable 
in our star due to non-LTE effects. Fe II lines, on the other hand, 
do not cover a large enough range of excitation potentials. We 
thus restricted ourselves to T^ff =6000, 6250, 6500 K and, us- 
ing the EW analyses of iron lines, deduced the following com- 
binations of logg/Vtur for these temperatures: 1.0/4.0, 1.5/4.0, 
2.0/4.0, with precision of 0.5 dex in log g (set by the model step 



size), and 1 km s"' for Viur- Our final best model of 6250/1.5/4.0 
is shown in Fig. [3] Linear fits to the Fell and Fel abundances 
plotted against EWs (and extrapolated to EW=0 for Fe I) give 
6.74 + 0.07 and 6.76 ± 0.06 for the iron abundance, respectively. 
The other two best-fit models are used for the estimation of the 
model-dependent errors on the abundances. 

3.3. Abundances 

The abundances of all elements deduced from the individual 
lines using our best model atmosphere are shown in Fig. |4] The 
elements are arranged in order of increasing condensation tem- 
perature to verify the presence of the depletion pattern reported 
in the literature for many disk sources. The average abundance 
per element and the error budget are shown in Fig. |5]and in Table 
[T] Given possible non-LTE effects, we assigned a quality flag to 
each element (each ionization state), that reflects the reliability 
of the final abundance estimate: flag 1 (best) - ions that do not 
show a trend with EW; flag 2 - neutrals for which at least three 
lines with EW<50 mA were available for averaging; flag 3 (least 
reliable) - all other cases. These flags are reflected in the circle 
sizes in Fig. |5] The errorbars are the RMS values, they result 
from the uncertainties in the EW measurement and the values of 
loggf-s. Table [T] also gives systematic errors due to the uncer- 
tainties in the model parameters. Singly-ionized species are only 
weakly sensitive to variations in Teff and Vtur compared to the 
neutral ones, but are more sensitive to log g. 

It is easy to see a moderate metal deficiency in BDh-46°442, 
at the level of -0.5 . . .-0.7 dex, compared to the Sun. Apart from 
the slightly more abundant C and S, however, there is no obvious 
trend with T^ond- The a-process elements (Si, S, Mg, Ti) show a 
0.1-0.3 dex enhancement over the iron group, while s-process 
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Fig. 3. Dependence of iron abundance in BD+46°442 on EWs 
and the lower excitation potentials of individual lines for our pre- 
ferred model atmosphere Teff/ logg/Vtm- - 6250/1.5/4.0. Teff 
was constrained from fits to hydrogen lines, Vtur from the re- 
quirement of no trend of individual Fe II abundances with EWs, 
and log g from the agreement between iron abundance inferred 
from Fell and weak Fel lines (strong Fel lines being more 
prone to non-LTE effects). The lower panel is less conclusive 
due to the non-uniform coverage of the excitation potentials by 
the Fe II lines. The dashed line indicates the iron abundance of 
the Sun. 



elements (Y, Zr, Ba, Nd) are not particularly enriched. The high 
abundance of N is likely a result of an evolutionary enrichment 
during the first and the second dredge-ups, but can also be af- 
fected by the non-LTE effects (Takeda 1992; Przvbilla & Butler 
12001 ). Overall this is a typical composition of a thick disk star 
without chemical depletion dReddv et al„.2006.) . 



4. Spectral energy distribution 

Fig. |6]compares photometric observations of BD-i-46°442 to the 
reddened photospheric Kurucz model with stellar parameters as 
derived in Sec. 13.21 The total reddening of E{B -V) = 0.18 was 
obtained by minimizing the difference between the adopted red- 
dened model and the observed optical and near-infrared fluxes. 
The presence of the excess emission beyond 2 /im is obvious. 
When we selected BD-i-46°442 for our sample, we could only 
rely on the IRAS detection of the excess (of which the 100 fim 
flux is still an upper limit), and now it is confirmed w ith the two 
other space missions, AKARI fehiharaetal1[20Toh and WISE 
(Wright et al. 2010). The SED, which is clearly not double- 
peaked and only slightly reddened, argues against a simple de- 
tached shell distribution for the circumstellar material. Instead, 
the SED of this kin d is characteristic of a star surrounded b y a 
passive disk (e.g, in'de Ru vter et al.ll2006l: lOielen et al.ll2007l) . 

We modeled the SED of BD-h46°442 wifli MCMax 
(iMin et alJ2009h . a 2D Monte Carlo radiative transfer disk code. 
This code computes the temperature structure and density of the 
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Fig. 4. Chemical composition of BD-i-46°442 as inferred from 
the EW measurements of individual lines in the 2009/08/17 
spectrum with the preferred model atmosphere Teff/ log §/ Vtur = 
6250/1.5/4.0. Colour-coded abundances are given relative to the 
Sun. Elements are arranged in the order of increasing condensa- 
tion temperature. 



disk. The vertical scale height of the disk is computed by an it- 
eration process, demanding vertical hydrostatic equilibrium. For 
the input stellar model we use the parameters as discussed in 
Section [3] and adopt a total system mass = 2Mq and a typ- 
ical post-AGB luminosity of = 2000 (corresponding to a 
distance for the system of 3.3 kpc). Since we lack information on 
the exact mineralogy of the dust in the disk, we use a mixture of 
astronomical silicates and metallic iron, with grain sizes ranging 
from 0.1 to 100 im\. The adopted gas-to-dust ratio is kept fixed 
at 100. The geometry of the disk is determined by the inner and 
outer radii Ri^ and /?out, the total dust mass Mjust, and the power 
law of the surface-density distribution of the dust E ~ R P. 

To keep the numbers of free parameters for the model low, 
we keep the outer radius fixed at 7?out = 500 AU, and use a value 
0.5 < p < 2.0 for the surface-density distribution, as expected in 
a disk environment. To account for the low observed total red- 
dening, the inclination is kept fixed at 15°. This modeling is quite 
degenerate and equally well fitting models with slightly differ- 
ent sizes, total masses and surface-density distributions can be 
found. Complementary measurements, such as interferometric 
observations and infrared spectroscopy, would be invaluable to 
constrain further the disk geometry and mineralogy. 

In Fig. |6] we show a disk model which gives good agree- 
ment with the photometric data. The disk parameters for this 
model ai-e: - 7 AU, Rom - 500 AU, p - 1.6, and Mdust = 
2 X 10"^ Mq. For the dust composition we use a combination of 
amorphous olivine and metallic iron, in a ratio of ~ 96% /4%. 
The AKARI and IRAS fluxes past 30 yum can not be fit by this 
model because that part of the SED requires even larger dust 
grains than adopted. These large grains are likely not uniformly 
mixed in the disk, but must have settled towards the mid-plane. 
Modeling such process is beyond the scope of our investigation. 

BD-i-46°442 thus belongs to a distinct group of evolved 
stars surrounded by dusty disks (Gielenet al. 2011). To verify 
whether it is a binary system, and whether the disk is circum- 
stellar or circumbinary, we analyze the radial velocities. 



5 



Gorlova et al: A new evolved binary with a disk 

Table 1. Chemical composition of BD+46°442. 
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Fig. 5. Abundances of BD-i-46°442 averaged from several lines 
of a given element, with error-bars representing the line-to-line 
RMS scatter. The circle size designates the reliability of the lines 
against the non-LTE effects (the bigger the more reliable, see 
text). No significant depletion pattern is observed for this star; 
week depletion is tentatively present, but better estimates are 
needed for the light elements to confirm it. 



5. RV curve 

We measured RVs using a cross-correlation method with a mask 
representing a G2 star The mask consists of ~1100 lines. Cross- 
correlation functions (CCFs) extracted from each order were av- 
eraged, and the resultant profile was fitted with a Gaussian to de- 
termine the RV. Some asymmetries in the CCF were noted that 
appear to correlate with the orbital phase (see discussion in the 
following section). They, however, were never large enough to 
affect the location of the main peak. The measured velocities are 
given in Table lA.ll and plotted in Fig.|7l Errors in the RV consist 
of two quadratically added terms: the error of the Gaussian fit 
and the 0.2 km s ' uncertainty due to the drift of the wavelength 
calibration zero-point. 

As one can see from Fig.|7] the three seasons of observations 
reveal a regular variation in the RV with a period ~100 days and 
peak-to-peak amplitude of ~40 km s"'. Integrating over half a 
period under the RV curve, we obtain a displacement of 130 Rq, 
which is more than the radius of a post-AGB star. The radial 
velocity variations cannot be due to pulsations. 

A long period with a nearly-sinusoidal shape of the RV 
curve, near-constant temperature, and the large displacement of 
the photosphere speak strongly against pulsations. Rather, we 
must be observing an orbital motion around an invisible com- 
panion. BD+46°442 is thus an SBl system. Indeed, we obtained 
a very good Keplerian solution with only one point (spectrum 
23) out of 60 deviating by more than 3cr (Fig. [8]). The derived 
orbital parameters are given in Table |2] The errors on the param- 
eters were estimated via the Monte Carlo method, using 1000 
simulated RV sets. For a given date a value for the RV was drawn 
from a Gaussian distribution with a mean and cr equal to the ob- 
served RV value and its error, respectively. For every simulated 
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Fig. 6. The observed spectral energy distribution of BD+46°442 
plotted against the reddened Kurucz's model (dashed line) and 
our star+disk model (solid line). The flux up to 30 fim can be 
reproduced with lO"*" Mq of dust in grains smaller than 100 //m 
in size that reside in a ring that stretches from 7 AU to ~500 AU. 



6 



Gorlova et al: A new evolved binary with a disk 




Fig. 7. HERMES radial velocities of BD+46°442 plotted against 
Keplerian orbit with parameters listed in Table |2] 



Table 2. Orbital elements of BD+46°442 
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Notes. Orbital parameters with their uncertainties and two statistical 
parameters indicating the goodness of the RV tit (reduced chi-squared 
and the coefficient of determination). 



dataset the orbital parameters were obtained, and their scatter 
was used as an error estimate on the parameters obtained from 
the real data. An example application of this procedure for the 
orbital period and eccentricity is given in Fig. |9] As one can see, 
albeit small (0.08), the eccentricity is significantly different from 
zero. 

Our observations provide the first account as far as we know 
of the RV variations in BD+46°442, that we successfully repro- 
duce with a Keplerian orbit. BD+46°442 proves to be yet another 
evolved system with a disk that is a binary. The semi-major axis 
of the primary is at most 1.2 AU (for / - 15°). This is smaller 
than the silicate dust sublimation radius that would be at 2 AU 
(given Li, - 2000 Lq), and much smaller than the actual inner 
disk radius at 7 AU estimated from the fit to the NIR excess 
with metallic iron particles. This barely leaves space for a cicum- 
stellar dusty disk around any of the binary components, unless 
the gas is very optically thick to shield some radiation. We con- 
clude that the dusty disk producing the IR excess is circumbinary 
rather than circumstellar. 



Fig. 8. Upper panel: radial velocities of BDh-46°442 phased 
with a period of 140.77 days and overlaid on our Keplerian fit. 
The horizontal line marks the corresponding systemic velocity. 
The (p - was defined to coiTespond to the maximum RV (i.e. 
near one of the elongations). Bottom panel: residuals from the fit 
with ±3cr margins. 

6. Line profiles 

6.1. Ha 

In Fig jA.ll we show that Ha alternates between two major pro- 
files: a double-peak emission with peaks separated by ~100 
km s"' and wings stretching to +200 km s"', and a P Cyg-like 
profile. This is better illustrated in Fig {10] (medium left panel) 
where we depict Ha after subtraction of the model photospheric 
profile. The resulting profile shows both emission and absorp- 
tion components that vary in strength with the orbital phase. 
The double-peak emission is observed most of the time, but is 
strongest and most symmetric between = 0.7 - 0.8, corre- 
sponding to the giant's inferior conjunction (when it is between 
us and the companion, Fig.fTTTi. When the giant starts to retreat, 
the emission diminishes, while the central absorption blue-shifts 
and broadens, quickly turning into a P Cyg profile. The expan- 
sion velocities reach 350 km s"' at(p- 0.3 - 0.4 when the com- 
panion is in front. Higher members of the Balmer series develop 
broadened cores at this phase, but the effect is much less pro- 
nounced than in Ha, particularly in the wings (Fig. lA.lt . which 
justifies their use for estimating Tgff/ logg . The Ha profile as a 
whole does not seem to shift in RV, but is not centered on the sys- 
temic velocity (Vsyst) either, rather, it is permanently blue-shifted 
by ~14km s"'. 

6.2. Non-photospheric components of metal lines 

In the other panels of Fig[TO](and in Fig. lA.2l l we show the CCF 
and a representative set of metal line profiles. As expected, the 
CCF follows the behaviour of weak to medium-strong lines, like 
Call 8912 A, that constitute the majority of unblended lines in 
the CCF list. These lines are "well behaved", being symmetric 
and hardly variable, except for the RV that reflects the orbital 
motion of the giant. The Fe 11(42) triplet, on the other hand, 
presents an example of strong lines (EW>300 mA) that show 
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Fig. 10. Dynamic spectra of selected lines and the CCF as a function of the orbital phase (time runs from top to bottom). Colours 
designate continuum-normalized fluxes (except for Ha), with black corresponding to the strongest absorption. Fluxes in the missing 
phases were obtained by a linear interpolation between the nearest observed phases. Dotted vertical line marks the systemic velocity, 
while the horizontal dashes on the left - the observed phases. One orbital period is shown twice to guide the eye. Ha is plotted after 
subtraction of the photospheric model spectrum; solid and dashed lines represent the RV curves (original and +50 km s ' off'set) of 
the giant and of the putative companion to illustrate the lack of significant motion in the emission component. Only data from the 
first two seasons of observations are shown (spectra 1-25, obtained between July 2009 - January 201 1); later observations confirm 
the general behaviour with phase, but the intensities of the circumstellar features appear to vary slightly from cycle to cycle. 
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Fig. 9. Distribution of the orbital period and eccentricity val- 
ues obtained from 1000 simulated RV data sets. The solid cross 
marks the solution based on the observed RVs (see Table |2]i, 
while the dotted cross and the upper right panel show an alterna- 
tive solution with a slightly better derived from the simulated 
RV data sets. 



an additional narrow component near Vsyst- It is best seen when 
the giant retreats from us ((f) = 0.9 - 1.1), while in the opposite 
phase it appears less prominent, being either weaker or blue- 
shifted (by ~10 km s"' from Vsyst) and overlapping with the pho- 
tospheric component. Other similar cases involve strong lines of 
Ball, Call, Fel, Fell, Mgl, and Till, and the resonance lines 
of K I. Exceptions include the prominent O I triplets at 7775 and 
8446 A, and the Call doublet at 8915 A, that still show nor- 
mal, single-deep, profiles. These transitions, however, require 
much higher excitation Ofiow =7-9 eV) and form deeper in 
the atmosphere. The strengthening of the central absorption be- 
tween <p = 0.8-1.25 and the disappearance or blue-shift between 
= 0.3-0.7 resemble the behaviour of the central absorption in 
Ha. Ha absorption, however, is more blue-shifted, particularly 
in the P Cyg phase. 

Nal D lines also show a central absorption component, but 
also an additional, narrower one 13 km s ' to the blue from Vsyst, 
which coincides with the centroid of Ha. Both components are 
probably of circumstellar origin, but it is difficult to study them 
due to blending with the photospheric lines. To the red, there are 
three supposedly ISM components, at RV= -23, -13, and -1 
km s that are stable in RV and depth. They are als o observed 
in the resonance KI lines. Using a code provided in lReid et alJ 
(l2009l) . we calculated kinematic distances to the foreground 
clouds responsible for these components. The largest distance 
is obtained for the most blue-shifted component (Vheiio = -23 
km s^', V['s'^ = -18 km s^'), which is 1.25 + 0.5 kpc. Taking 
this as a minimum distance to BD+46°442 itself, and correct- 
ing for the reddening Ay = 0.56, we obtain a lower limit on the 
stellar luminosity: My = -1.54^Qgp Clearly, the star is brighter 
than normal EG giants of the luminosity class III (My = 1 ... 2), 



but is consistent w ith being an F type post-AGB star, which have 
Mv = -2 . . . - 3 (iBondll 19971: iGinestet & Carquillatll2002 V On 
the other hand, the two components at Vheiio ~ -100 km s"' 
would indicate an lab supergiant (My = -6 ... - 7) at more than 
2.5 kpc above the Galactic plane, which is very unlikely. Thus, 
the two components near Vsyst are most likely circumstellar. 

In addition to the absorption components, as can be seen 
from Eig.[TO] near =0.75 corresponding to the giant's inferior 
conjunction, weak emission appears in the wings of the sodium 
lines, echoing Ha. We observe only one other system of emis- 
sion lines beside Ha and Nal D: the weak Eel(12) multiplet at 
8047.6 and 8075.1 A. Similar to the central absorption of metal 
lines, the emission is centered on the systemic velocity and does 
not participate in the orbital motion (Eig.fTOli. 

6.3. Secondary spectrum ? 

The metal lines noticeably strengthen near the giant's inferior 
conjunction (0 = 0.75). Zooming in on the wings of the CCE 
(see the right-bottom panel of Eig. [TOl i. however, reveals also 
faint absorption spectrum moving in anti-phase with the gi- 
ant's spectrum. Therefore, we may be observing a companion 
spectrum instead. We remodeled the CCE profiles with two 
Gaussians, where possible, and obtained a semi-amplitude of 
15+5 km s ' for the secondary spectrum. Dividing by the 23.7 
km s"' of the primary, we obtain a mass ratio of 0.6. Given this 
mass ratio and the mass function of 0.2, we obtain the individ- 
ual masses of 0.9-0.3 Mq for the giant and 1.4-0.5 for the 
companion for a range of inclinations between 45-90°. These 
masses are consistent with a low-mass pair where the originally 
more massive component has evolved into a post-AGB, while 
the companion is still on the MS or in the sub -giant stage. Usin g 
expressions for the Roche-lobe radius from 'Paczvriskil (Il97lh . 
that depend only on the mass ratio, we obtain for the giant 
Ri ~ 0.37 X (fli + fl2) = 57 - 39Re. The radius of BD+46°442 
matches this range well if it has a typical post-AGB luminos- 
ity of 2000-3000 Lq. Therefore, we are observing an interacting 
system where the post-AGB star is (close to) filling its Roch-lobe 
and as such feeding a more compact companion. 

We want to stress, however, that these are very preliminary 
estimates. The depth of the secondary component in the CCE is 
about 1% (compared to 10% of the giant's spectrum), which is 
far below the typical S/N ratio of our spectra. The presence of 
this component thus can not be verified from individual lines. 
Furthermore, the spectral type of the companion can not be sig- 
nificantly different from that of the primary, otherwise one would 
not detect it with the G2 mask. The companion then should 
be only ten times fainter than the giant primary, which is too 
bright were it a MS star. And if both components are giants, 
they would have to have nearly identical initial masses due to the 
short longevity of the giant stage. It is possible instead that we 
are seeing a pseudo-photosphere of an accretion disk around the 
companion rather than the companion's spectrum itself. Higher 
S/N spectra and a multi-band photometry are needed to confirm 
the contribution of the companion in the system's flux. 

7. Discussion 

The discovery of the RV variations in BD-i-46°442 provides a 
strong support for the disk-binary connection in post-AGB stars. 
With a period of 140.77 ± 0.02 it is one of the shortest known 
post-AGB binaries. A small but significant eccentricity indicates 
that the system is not completely circularized. In the following 
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discussion we will show how the observed spectral features in 
BD+46°442 indicate an ongoing mass transfer from the giant to 
the companion via the Roche-lobe overflow. 

7.1. Ha variability in interacting binaries 
The Ha pro files discussed here are not uncommon among post- 



AGB stars (IWaters et al.lll993bt iP ollard e t al.lll997t iMaas et alJ 
120051; ISanchez Contreras et al.T l^OOS), but they have not been 
systematically studied for variability in function of orbital phase. 
We are aware of four other disk systems where a strong blue- 
shifted absorption was noted to develop during the giant's su- 
perior conjunction: H R 4049 with Pgrb = 430*^, a sin / - 0.6 
AU, > 60°, e = 0.3 (IWaelkens etanil99Ul>^n Winckel et al.1 



119951; lBakker etal.1 119981) ; HP 44179 (the central star of the 
Red Rectangle nebula) wifli P^.b = 318'', a sin; ^ 5 AU , 
effective" = 35°, 6 = 0.4 (IWaclkens et al.lll996l;IWitt etani2009h ; 
IRAS 8544-4431 wifli ^ 508**, a sin;' ^ 0.4 AU, ; ~ 60° 
e ^ .2 dMaas et alj|2003l; iDeroo et"ani2007l: Ivan Winckel et all 
12009^; IRAS 19135+3937 with P^^t = 127^ a sin/ = 0.2 AU, 
e ~ 0.3 (G orlova et al.ll2012l) . The fact that BD-h46°442 has a 
very small eccentricity (e = 0.083 + 0.002) indicates that the 
cause of the spectral variations most likely lies in the varying 
line of sight towards the components, rather than in the physi- 
cal changes during the orbital motion, such as an increased mass 
loss at periastron. 

The most common interpretations of the double-peak Ha 
emission in stellar spectra are: an inclined Keplerian gaseous 
disk (e.g. in Be and T Tau stars), the giant's atmosphere 
that is irradiated by a hot companion (Cataclysmic Variables, 
Symbiotics), propagating shock in pulsating stars (Miras, RV 
Tau stars), and the chromosphere. There is no indication in our 
spectra of a hot companion (no helium lines or nebular emis- 
sion), pulsations are excluded based on the smooth orbital RV 
curve, and the existence of chromospheric activity has yet to be 
established in post-AGB stars. We therefore concentrate on the 
disk hypothesis. The emission could originate in the gaseous ex- 
tension of the circumbinary dusty disk towards the center of the 
system or in a circumstellar disk around one of the components. 
To differentiate between these possibilities we searched for sim- 
ilarities with other types of interacting binaries. There, a disk 
sometimes forms around a more compact companion as a result 
of the Roche-lobe overflow or a wind accretion from a larger 
companion. 

In Algols, when the binary separation is large, the accretion 
stream from a cool giant curves and settles into a disk around a 
more massive hot MS primary (Richards & Albright 1999). The 
simulations of iMiller et al.l(l2007i) show that the dominant source 
of Ha emission in Algols is the disk, while the stream provides 
a much smaller contribution. A similar situation could be oc- 
curring in BD+46°442, only the accreting companion is much 
fainter than the giant (Fig. fTTT i. The problem with this interpre- 
tation of the Ha emission is the lack of convincing motion in 
the anti-phase with the absorption spectrum of the primary (Fig. 
[Tol l. The P Cyg-like profile, on the other hand, is a signature 
of an outflow. It is very rare to find such broad absorption fea- 
ture (with velocities up to -300 km s~') in Algols, and it may 
not be observed in every cycle (IPeterslll989[ iMiller et al.ll2007h . 
Usually Algols display classical and inverse P Cyg-like profiles 
at a smaller range of velocities just before and after giant's in- 
ferior conjunction, which is explained by the disk eclipse. This 
is clearly not the case in BD+46°442 where the wind profile is 
strongest near the giant's superior conjunction. The explanation 
of the P Cyg profile, therefore, should be searched elsewhere. 



W Ser systems are thought to be in the stage preceding 
Algols, when the transfer rate is very high and the accretion disk 
completely obscures the a ccretor In at least on e such system, 
V Sgr with Poi-b = 137'* (iNetolickvet al.ll2009l) . a P Cyg-like 
profile is obse rved near the giant's superior conjunction, which 
iNariail (Il967l) explained by a "coronal stream" from the giant 
to the companion. At this orbital phase the stream is directed 
towards us and at the same time is projected against the giant 
(Fig.fTTTi. creating a blue-shifted absorption. Furthermore, in the 
simulations of interacting binaries the accretion stream is often 
found to be curved in the direction of the giant's motion, which 
could explain why the largest speeds in BD+46°442 are reached 
at a slightly later phase after the conjunction (compare spectra 6 
and 7 on Fig. lA.lb . Interestingly, v Sgr is considered to be a non- 
eclipsing type of W Ser stars. In classical, edge-on, systems the 
blue-shifted absorption is not observed, perhaps because the ac- 
cretion stream, when projected against the donor, is hidden from 
our view by the accretion disk. 

Symbiotic stars are another relevant class of objects to con- 
sider. There double-peaked emission profiles are quite common 
( Van Winckel et al. 1993), and are being interpreted as either 
due to an accretion disk around a WD (Robinson et al. |1994|: 
IWn Eck & Jorissenll2002l) or a disk wind (Skopal.2006.) . There 
are even reports that the blue peak gets suppressed at the orbital 
phase corresponding to the occurrence of the P Cyg-like pro- 
file in BD+46°442 (Robinson et al. 1994; Murset et al. 2000l), 
but the true broad absorption in symbiotics is only observed in 
outbursts. 



7.2. Jet iauncliing in the accretion disk of the companion 

Summarizing, the variations of Ha in BD-i-46°442 are complex, 
but can be decomposed into several time-dependent properties 
for which we need an explanation. 1) The double-peak emission 
is rather stable in velocity along the orbit, which means that it is 
connected to the system itself and not to one of the components. 
The emission, however, varies in strength, being strongest dur- 
ing the giant's inferior conjunction. 2) The P Cyg-like absorption 
component is observed at very high velocity, but only during a 
very short orbital phase interval, when the giant is at the superior 
conjunction (see FigfTOt. As we have seen in the previous sec- 
tion, some of these aspects can be explained by the interaction 
with the companion. 

It is useful to compare this behaviour to the line va r iabilit y 
as observed and discussed in HD 44179 bv I Witt et all ( l2009l) . 
IWitt et all (l200 9») attributes both emission and absorption in Ha 
to the same structure, which is a jet launched from the accre- 
tion disk of the companion. P Cyg-like profile forms when the 
blue-shifted lobe becomes projected against the face of the gi- 
ant, which is impossible for the red-shifted lobe and which is 
therefore always observed in emission. The jet launched in a pre- 
cessing accretion disk around the comp anion of HD 44179 may 
have created the Red Rectangle nebula (IVelazquez et al.ll201 ll) . 
This model can be adopted for BD-i-46°442 (Fig.E), with one 
caveat. Red Rectangle is a peculiar system, because the optical 
flux is dominated by scattering as our line of sight is in the or- 
bital plane of the binary, allowing to see both jets at a range of 
inclinations. For BD-i-46°442, on the contrary, the smaller IR ex- 
cess indicates that we are observing the star directly, along the 
line-of-sight which is not in the orbital plane. The small line- 
of-sight reddening shows that our aspect angle towards the sys- 
tem is not through the puffed-up circumbinary disk. The blue ab- 
sorption as seen at conjunction is created by the blue-shifted jet 
which is projected on the giant only in this specific orbital phase. 
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Continuum photons coming from the giant will be scattered out- 
side the line-of-sight by Ha resonant scattering on high velocity 
hydrogen atoms in the jet. The Doppler shift of 300 km s ' cor- 
responds to the line-of-sight velocity component and needs to be 
deprojected to obtain the outflow velocity of the jet. As we do 
not constrain the jet opening angle, nor the angle between the jet 
axis of symmetry and the orbital plane, the deprojected veloc- 
ity is difficult to quantify. The static double-peak emission still 
requires a different origin, as discussed in the previous section. 

The existence of jets has been strongly advocate d based on 
other i ndependent argu ments in e.g., a /? Lyr system (iHarmaned 
120021: lAk et al.l |2007'), a proto-type of the massive Roche- 
overflow binaries (of which W Ser is a sub-class). Jets have been 
invoked to explain bipolar planetary and p roto-planetary neb- 
ulae (PPNe), including the Red Rectangle dCohen et alJl2004 
IVelazquez et al. 2011). Sahai et al. (2002) detected a proper mo- 
tion of jets in the PPNe Henize 2-90 with velocity of a few hun- 
dred km s which is consistent with the Keplerian velocity at 
the surface of a low- to intermediate-mass MS star. Finally, jets 
have been also mentioned in symbiotics, in regard to symmet- 
ric bu mps at +1000 km s"' appearing in outbursts (Skopal et al] 
12009 ). though they have never been invoked to explain the 
double-peaked profile of the quiescent state. It is difficult to in- 
clude jets in the hydrodynamic calculations due to the supersonic 
velocities, or in the radiative transfer calculations. Nevertheless, 
the possibility of jet contribution int o Ha and HP formation 
was demonstrated in several studies (iBudai &Richardsl 12004 
lArrieta et al1l2005l) . 

7.3. Gas stream 

Further phenomenological similarity between BD-i-46°442 and 
the mass-exchanging binaries can be seen in the non- 
photospheric "shell" components of metal lines. The narrow 
features, mostly in absorption, discussed in BD+46°442 (such 
as the Fe II (42) and the near-I R Ca II triplets), have been first 
noted in fast rotating hot stars (ISlette bak 1986). They origmate 
from the meta-stable levels at 2-4eV above the ground state, 
and therefore must form in a low-density, but hot (several kK) 
circumstellar medium. In cooler stars, shell components have 
been mostly studied in the resonance Call H&K and Nal D 
lines, but these transitions are prone to interstellar contamina- 
tion. Because of their regular behaviour with the orbital phase 
in BD+46°442, we again consider them to originate in the bi- 
nary flows, rather than in some clumps in the circumbinary disk 
(an expla nation proposed for the protoplanetary disks, see e.g. 
iMora et al. (2004)). The behaviour of these features appears to 
be very diverse from system to system, so it is not surprising 
that they were ascribed to a whole range of accretion-related 
struc tures: disk, disk-wind collision shocks, jets, streams, etc . 
(e.g. I Andersen et al.' '1988'; 'Plavec' '1988'; 'Weiland et al.' '1995"; 
[Miller et al. 2007; Sudar et al. 201 1; Mennickent & Smith 2010; 
iHarm anec 2002; Ouiro ga et al.ll2002l) . 

We found a particularly good agreement between the be- 
haviour of the metal absorptions in BD+46°442 with He I lines 
in a W Ser-type star RY Set (Porb = 1 1'', orbital separation ~0.2 
AU), as described in Grundstromet al. (2007). To explain these 
features, one needs gas projected on the giant over at least half 
of the period and moving with a slower velocity, for example a 
circumbinary Keplerian disk. Furthermore, this ring must have 
an asymmetric density distribution to explain unequal strength 
of the absorption in elongations. Grundstromet al. (2007) pro- 
posed that this structure could in fact be a plume on the trail- 
ing side of the giant, an outflow from the L2 point that ex- 




Fig. 11. Schematic representation (not to scale) of the views of 
BD+46°442 at different orbital phases. The dusty circumbinary 
disk has been omitted to better illustrate gas flows in the vicin- 
ity of the components (the inner disk wall would be outside the 
boundaries of the figure anyway). The inclination, both of the 
orbit and of the likely precessing jet, remains to be constrained. 
The arrow shows the direction of the giant's motion. The figure 
illustrates how an accretion jet around the companion can ex- 
plain the transient P Cyg-like profile of Ha, while a gas stream 
trailing behind the giant can explain the narrow absorption com- 
ponent of metal lines. 

tends beyond the orbital plane and warps around the system 
(Fig. [TTT i. It explains why the central absorption is clearly visi- 
ble when the giant recedes from us, but is much less pronounced 
in the opposite elongation. This outflowing arm is consistently 
reproduc ed in the hydrodynamica l simulations of interactive 
binari e s (Ide Va l-Borro et al. 20091 iMohamed & Podsiadlowskil 
1201 Ih . iGrundstrom et al.l (I2007i) further speculate that the out- 
flows through the L2 and L3 points could explain the formation 
of the large nebula (R~2000 AU) surrounding RY Set. This sce- 
nario brings an intriguing possibility of the current formation of 
the circumbi nary disk aro und BD-h46°442. Finally, in a recent 
studv. lThomas et alJ (|201 II) reports on the identical behaviour of 
some metal lines, includ ing Fe II 4924 A , in HD 44179. Unlike 
IGrundstrom et al.l (l2007l) . they interpret it in terms of the photo- 
spheric line asymmetries due to matter outflowing from the giant 
near periastron and falling back near apastron. This model, how- 
ever, is questionable for BD+46°442 due its low eccentricity. 
Obviously, observations of more systems with different eccen- 
tricities and orbit orientation are needed to test these hypotheses. 

Another manifestation of the circumstellar gas around 
BD+46°442 are two weak emission lines of the Fel(12) multi- 
plet in the far red region of the optical spectrum. These emis- 
sion lines have been reported at least in two other disk post- 
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AGB sources, 89 Her (a binai-y with foA = 288'', a sin / = 0.1 
AU, i - 15°, e - 0.2, and a resolved bipolar CO outflow, see 
IWaters et ai] |1993a; Buiarrabal et al. 2007) and a candidate RV 
Tau star QY Sge (Kameswa ra Rao et alJl2002h . In 89 Her these 
lines are narrow (FWHM<10 km s"') and single-peaked, while 
in QY Ser and BD+46°442 fliey are much broader (FWHM~40 
km s"') and perhaps even double-peaked. Unfortunately, the 
latter is difficult to measure accurately for BD-h46°442 due 
to con tamination with telluric absorption. iKameswara Rao et alj 
(l2002h proposed that this emission originates in the bipolar wind 
from the system, in particular, from the gas within the inner cav- 
ity of the circumbinary disk that re-emits starlight scattered by 
the disk walls. Assuming a total stellar mass for BD+46°442 of 
1 .5 Mq and a circumbinary gas disk with Keplerian rotation, the 
HWHM=20 km s ' indicates a distance of 3.3 AU. Given that 
the inner radius of the dusty disk is at 7 AU, as deduced from 
the SED fitting, this gas ca n indeed be located within the inner 
disk hole. .Kameswara Rao et alj (l2002h also explained broader 
emissions in Ha and the Na I D lines as originating in the same 
wind, only further out from the disk plane where the wind has 
been accelerated. This scenario would explain the lack of the 
RV variations in the emission lines of BD-i-46°442, including 
Ha. The remaining differences between 89 Her and QY Sge 
with BD-i-46°442 could be due to different inclination angles. 
To constrain the latter for BD-i-46°442, an accurate light-curve 
is needed. 

The next obvious step to confirm the association of all these 
circumstellar features with binarity, would be a systematic com- 
parison of disk sources to the sources with spheri cal shells (as 
inferred from the SED, e.g. I van Aarle et all (1201 lb ), as the latter 
presumably hai^bor single stars. 

8. Summary 

We obtained 60 echelle spectra over a period of 3.5 years for 
BD-i-46°442, a poorly-studied high-galactic giant with a dusty 
disk. We derive the following photospheric parameters: T^ff - 
6250 + 250 K, logg = 1.5 ± 0.5, and an average metallicity 
[M/H] = -0.7 + 0.2, without a strong depletion pattern. The 
enhanced abundance of a elements is characteristic of the origi- 
nal composition of a thick-disk low-mass star, consistent with a 
post-AGB interpretation. The observed large amplitude of the 
RV variations and a lack of a strong variability in Teff argue 
against pulsations. We therefore attribute the RV variations in 
BD+46°442 to binarity (with a tentative detection of a compan- 
ion spectrum in the CCF). We find the following values for the 
orbital parameters: an orbital period 140.8 days, an eccentric- 
ity 0.08, and separation <1 AU. This adds to the several dozen 
post-AGB disk binaries with known orbital parameters. The or- 
bital period falls on the short side of the period distribution of 
the other orbits. 

Time-resolved spectroscopy allowed to detect gas streams in 
BD+46°442, indicating that it is an interactive binary. Ha (and 
to a smaller extent higher Balmer lines and the Na 1 D) are found 
to alternate between double-peaked emission, which is charac- 
teristic of a Keplerian disk, and a P Cyg-like profile, which is 
characteristic of an outflow. The blue absorption reaches values 
up to 300 km s ' and develops only around the giant's superior 
conjunction. We suggest it is due to a jet that originates in the 
accretion disk around the companion. In addition, strong metal 
lines display central or slightly blue-shifted absorption inside the 
photospheric profiles. Similar features are also observed in some 
massive interactive binaries and likely in several other post-AGB 
disk sources (like in the famous Red Rectangle). 



The discovery of mass transfer and, possibly, accretion- 
driven jets in post-AGB binaries opens a new avenue to ex- 
plore some of the long-standing puzzles of these objects, such 
as a variable mass-loss rate, a long life-tim e of the circumbinar' 
disk s, a lack of exte nded dusty nebulae tSiodmiak et alj [200 
Lag adec et alj 201 1[). and the b ipolar structure of the gaseous 
outflows (Buian'ab al et al.ll2007l) . By uncovering more systems 
like this in the course of our spectroscopic survey, aided by com- 
plimentary observations, we hope to pin-point the exact geo- 
metrical structures responsible for the various non-photo spheric 
features. This will enable to construct realistic hydrodynamical 
models, providing a more solid basis for the theory of binary 
evolution. 
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Appendix A: Representative line profiles 

Fig. lA. II shows Ha and three other Balmer lines in BD+46°442 
ordered according to the orbital phase, where each of the five or- 
bital cycles (covering the first 1.5 years of observations) is desig- 
nated with a different colour A corelation with the orbital phase 
is obvious. The central absorption in Ha is not always aligned 
with the photospheric velocity, which means that the profile can 
not be represented by a simple superposition of a broad emission 
with a photospheric absorption. Near giant's inferior conjunction 
((f) - 0.75) Ha exhibits a double-peak emission, while near su- 
perior conjunction (0 = 0.25) the blue peak is replaced with 
an extended blue absorption, as can also bee seen in the higher 
Balmer lines. 

Fig. IA.2I shows the behaviour of some metal lines and the 
CCF The CCF represents the behaviour of the majority of the 
photospheric lines, that are symmetric and move with the or- 
bital velocity. Low-excitation (xiow -0-3 eV), strong (EW >300 
mA) lines, on the other hand, show additional circumstellar 
and interstellar (in Nal D) components, particularly between 
(f> = 0.9-1.1. 
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Fig. A.l. Balmer lines as a function of the RV phase (0 = corresponds to the maximum redshift). Different colours denote different 
orbital cycles. A solid vertical line marks our systemic velocity of -98.9 km s"', while dotted lines mark 50 km s"' intervals from 
it. The arrows mark the photospheric velocity according to the CCF. 
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Fig. A.2. Same as in Fig lA.ll only for the CCF and the representative metal lines of different strengths: Call 4d-4f doublet, Fe II 
(42) triplet, and the Na I D doublet. In each panel different line styles (solid, dashed, dotted) designate different members of the 
same multiplet. 
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Table A.l. Radial velocities of BD+46°442 



N 


TT^ Acc r\f\r\ , 

JD 2,455,000+ 


RV 


(T 




(d) 


(km s ') 


(km s ') 


1 


31.71 


-89.62 


0.26 


2 


42.70 


-80.72 


0.30 


3 


60.70 


-75.28 


0.33 


4 


60.72 


-75.29 


0.33 


5 


85.64 


-92.73 


0.27 


6 


92.64 


-100.15 


0.28 


7 


100.62 


-109.71 


0.28 


8 


125.49 


-122.42 


0.32 


9 


131.53 


-120.76 


0.30 


10 


156.40 


-102.56 


0.26 


11 


159.55 


-100.60 


0.26 


12 


159.56 


-100.42 


0.27 


13 


199.33 


-75.91 


0.30 


14 


216.38 


-83.05 


0.29 


15 


224.36 


-90.37 


0.29 


16 


421.63 


-116.87 


0.30 


17 


425.67 


-112.94 


0.30 


18 


432.59 


-106.89 


0.27 


19 


468.66 


-78.39 


0.30 


20 


471.62 


-76.86 


0.32 


21 


499.54 


-83.24 


0.26 


22 


507.47 


-91.20 


0.27 


23 


553.48 


-123.07 


0.27 


24 


573.34 


-107.14 


0.25 


25 


581.34 


-100.30 


0.22 


26 


767.71 


-76.56 


0.30 


27 


769.71 


-77.36 


0.30 


28 


776.73 


-80.47 


0.29 


29 


783.69 


-85.43 


0.27 


30 


791.74 


-94.47 


0.25 


31 


797.73 


-101.45 


0.26 


32 


826.65 


-123.22 


0.28 


33 


835.57 


-121.05 


0.29 


34 


841.50 


-118.54 


0.28 


35 


842.59 


-117.29 


0.30 


36 


845.61 


-114.46 


0.29 


37 


848.41 


-112.60 


0.28 


38 


852.52 


-109.55 


0.25 


39 


857.64 


-104.90 


0.26 


40 


860.60 


-102.75 


0.25 


41 


866.54 


-98.31 


0.25 


42 


870.56 


-95.45 


0.25 


43 


872.61 


-92.87 


0.24 


44 


880.56 


-84.86 


0.26 


45 


882.46 


-83.74 


0.26 


46 


884.55 


-82.57 


0.27 


47 


886.49 


-81.61 


0.26 


48 


888.55 


-80.62 


0.25 


49 


890.52 


-79.87 


0.25 


50 


890.54 


-79.79 


0.25 


51 


903.38 


-74.61 


0.25 


52 


911.41 


-77.25 


0.28 


53 


934.41 


-96.74 


0.30 


54 


940.36 


-105.17 


0.31 


55 


943.44 


-108.88 


0.30 


56 


948.40 


-114.77 


0.28 


57 


953.36 


-118.43 


0.27 


58 


956.40 


-119.74 


0.26 


59 


957.48 


-120.08 


0.26 


60 


968.45 


-121.61 


0.33 



Notes. Radial velocities (barycentric correction included) measured 
with HERMES over the period July 2009 - February 2012. 



